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Lecture notes presented at Windsor NATO school on " Field Theory of Strongly Correlated Fermions 
and Bosons in Low-Dimensional Disordered Systems" (August 2001). The purpose of these lectures 
is to give a brief modern introduction to Keldysh non-equilibrium field theory and its classical 
analog - Doi-Peliti technique. The special emphasis is put on stressing the analogy between the 
two approaches. 



I. QUANTUM KELDYSH FORMALISM 

1. A hit of history. The name of the technique takes its origin in 1964 paper of L. V. Keldysh . Among the earlier 
closely related approaches, one should mention J. Schwinger 1^ and R. P. Feynman and F. L. Vernon g]. Classical 
counterparts of the Keldysh technique are extremely useful and interesting on their own. I shall mention Wyld 
diagrammatic technique [j7|, which is the basis of the modern theory of turbulence; Martin-Siggia-Rose-DeDominics 
method 1^ for stochastic systems; and Doi-Peliti technique ]9|JlO|] for reaction-diffusion systems and cellular automata. 
Despite the unique power and wide popularity of the technique its pedagogical presentations are rare. One may find 
some in books or in a review This lecture is not an alternative for the existing literature. Using the simplest 
possible example, I try to explain the general structure of the theory, its relation to other approaches, its potential 
for novel applications and open problems. Technicalities of the numerous practical applications may be found in the 
literature and left outside the scope of the lecture. I tried hard, however, to mention possible traps and silent points 
which are rarely exposed explicitly. 



2. Motivation. The Keldysh formulation of the many-body theory is useful in the following cases: 

• For treatment of systems not in the thermal equilibrium [Q . 

• For calculation of the full counting statistics of a quantum mechanical observable (as opposed to an average 
value or correlators) |]ll| , |l2| . 

• As an alternative to the replica and supersymmetry methods in the theory of disordered and glassy systems 



For equilibrium problems where Matsubara analytical continuation may prove to be sophisticated 



3. Closed time contour. The standard construction of the zero temperature (or equilibrium) many-body theory 
(see e.g. [|^D) involves the procedure of adiabatic switching interactions on at distant past and then switching it off 
at distant future. The crucial assumption here is that starting from the ground (or equilibrium) state of the system 
at t = — oo one arrives to the same state a,t t = +00 (acquiring some phase factor along the way). This is clearly not 
the case out of the equilibrium. Starting from some arbitrary distribution and then switching interactions on and off, 
one is going to find the system in some unpredictable state. The latter depends, in general, on the peculiarities of the 
switching procedure. The absence of knowledge about the finite state spoils completely the entire construction JTs],^. 
One would like thus to build a theory which avoids reference to the state aX t = -\-oo. We still need to know a finite 
state (since we calculate traces). The Schwinger suggestion is to make the finite state to be exactly the initial one. 
The central idea is to let a quantum system to evolve first in the forward direction in time and then to "rewind" its 
evolution back, playing the movie in the backward direction. One ends up thus with the need to construct a theory 
with the time evolution along the two-branch contour depicted on Fig. |l|. Then, no matter what was the state at 
t — +CXD, after the backward evolution it will return back to the known initial state. In this construction there is no 
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switching of interactions in the future. Both switchings on and off take place in the past: on - on the forward branch 
and off - on the backward. How to construct such a theory and how to use it - is the subject of this chapter. 
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FIG. 1^ The closed Keldysh time contour C. Dots on the forward and the backward branches of the contour denote the discretized 
time points. 



4. Field theory. To be specific, let me consider the simplest possible many-body system. It consists of bosons living 
in a single quantum state with the energy ojq: 

H ~ ujQa)a ; (1) 
here and a are bosonic creation and annihilation operators. Let us consider the "partition function" defined as 

Z = Tr{pof/c}/Tr{po}, (2) 

where Uc is the evolution operator along the closed contour C. If we assume (as we shall do for a while) that all 
the external fields are exactly the same on the forward and backward branches of the contour, then J7c = 1 and 
therefore Z ~ 1. In Eq. po = p{H) is some |^ density matrix operator defined sX t — — oo. In our example 
Tr{po} = [1 ~ ^('^o)]"^- An important point is that in general Trjpo} is an interaction- and disorder-mdependeni 
constant. Indeed both interactions and disorder are supposed to be switched on (and off) on the forward (backward) 
parts of the contour sometime after (before) t ^ —oo. I shall therefore frequently omit this constant or refer to it as 
being unity - it never causes a confusion. 

The next step is to divide the C contour into 2N + 1 time steps of the length 5t, such that ti — t2N+i = ^oo and 
tN+i = +00 as shown on Fig. |^. Following the standard route we obtain the coherent state functional integral, 
by introducing a resolution of unity at each time step. Taking the N oo limit we obtain for the partition function 

Z = I e^^{iS[m,m]} , (3) 



Tr{po} 



where the action is given by 



Sm,m = J dt^{t)V-'<t>{t), (4) 
c 

and 'D~^ — [idt — luq]. It is important to remember that this continuous notation is only a short way to represent the 
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2N X 2N matrix Vr^ 



^To accommodate the equilibrium initial density matrix, po — exp{— /9(i/ — f^N)}, it is sometimes suggested to add an 
imaginary (vertical) part to the Keldysh contour a,t t — —oo. To my opinion this procedure only obscures the structure of the 
theory. Indeed, the technique is not restricted in any way to the initial distribution being the equilibrium one. Therefore it is 
much preferred to keep the structure general and not commit oneself with unnecessary assumptions. 



2 



1 

-1 + StiLUo 




1 

-1 + StiujQ 



(-1 + Stiu;o)p{uJo) 




1 

-1 + Stiujo 



-1 - 6tiuJo 



1 

-1 — Jfici^o 



(5) 



The change of sign +6t — > —St takes place in the center of the matrix at the point where the time axis changes its 
direction. It is straightforward to evaluate the determinant of such matrix 



det[ - iV-'] = 1 + (-l)2^+V(^o)(-l + StiLJo)''{-l - StiuJof ^ 1 - p(c^o) 



(6) 



Recalling that p{uJo) < 1, one may convince himself that all eigenvalues of the matrix —iD^^ have positive real part. 
This fact secures convergence of the functional integral. Performing the integration, one not-surprisingly obtains for 
the partition function 



det"^ [ - i2?-i] 
Tr{po} 



(7) 



Seemingly absent in the continious notations - the upper right element of the matrix is extremely important to obtain 
the correct normalization. It is not less important to find the proper correlation functions. 

We divide next the bosonic field (f){t) onto the two components 4>f{t) and (t>b{t) which reside on the forward and 
the backward parts of the time contour correspondingly. The action may be now rewritten as S'[(/)/] — >5'[(/)f,], with the 
time integration running along the real axis (the minus sign comes from reversing direction of the time integration on 
the backward part of the contour). 



5. Correlation Junctions. Now comes the most unpleasant part: to invert the matrix and find the correlation (or 
Green) functions. We may save ourselves a lot of efforts by recalling that the functional integral always produces the 
time-ordered Green functions |l| ] and that the time arguments on the backward branch are always after those on 
the forward. This way one findsT 

- i{<l)f{t)4>f{t')) = V^{t,t') = 9{t - t')V>{t,t') +0{t' - t)V<{t,t') ; 
-z{cj,f{t)Mt'))='D<{t,t') = (-*ne'*'^«(*-*')) ; 

-j(0b(t)^/(O) =2?>(t,0 = (-*("+ 1)6-^""^*"*'^' 
-i{(l)b{t)4>f{t')) = V^{t,t') = 9{t' - t)V>{t,t') + 0{t - t')V<{t,t') , 



(8) 



where the angular brackets are understood as the functional integration with the action given above. The structure 
is general, but 2?^ = —i{\a{t)a'^ {t')\) and — —i{\a^t')a{t)\) are given for orientation in brackets for our simple 
example (n = p{uJo)/{l — p{oJo)) is the bosonic occupation number). 

It may be surprising to see the non-zero off-diagonal correlators {4'f4'b) and {(j)h4>f), while the action seems to be 
perfectly diagonal in the f — b space. The trick is exactly in the presence of the upper-right element of the discretized 
matrix, hidden in the continuous representation. Notice, that the presence of the off-diagonal elements does not 
contradict to the continuous notations. Indeed, [idt — wojl^^'^ — 0, while [ii9t — ^q\D^-^ — ±(5(i — t'). Therefore in 



^The more familiar definition of the Green functions is D = ~i{(f>(j)) / Z , since in our case Z = 1 - we are doing OK. 
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the obvious 2x2 matrix notations T>~^'D = 1, as it should be. The point is that the [idt — operator is not 

uniquely defined. The boundary conditions must be specified and the upper right element does exactly this. 
Obviously not all four Green functions defined above are independent. Indeed, an inspection shows that 

+ +V< . (9) 

One would therefore like to perform a linear transformation of the fields to benefit explicitly from this relation. This 
is achieved by the Keldysh rotation. 



5. Keldysh rotation. Let us define the new fields as 

Mt) = \{<l>f{t) + Mt))- ci,,{t) = ^{<i>f{t)-Mt)) (10) 

with the analogous transformation for the conjugated fields. The integration measure changes by unessential multi- 
plicative constant. The subscripts cl and q stand for classical and quantum components of the fields correspondingly. 
The rational for such notations will become clear shortly. First, a little algebra shows that |^ 



where a, (3 = cl, q. Here indexes i?. A, K stay for retarded, advanced and Keldysh components of the Green function. 
These three are the main objects of the Keldysh technique. They are defined as 







-V< 


= 9{t- 


-t'){V> -v< 






-V> 


= e{t' 


-t){V< -v> 




= V> 


+ V< 







(12) 

Notice that the retarded (advanced) component is lower (upper) triangular matrix in the time space. Moreover, our 
regularization of the functional integral guarantees that these functions never appear at exactly coinciding arguments, 
removing thus uncertainty about the value 9{t — t) = 9{Q). Fields (j) always appear at one time step 5t after (j) fields 
on the Keldysh contour. As a result, 0/(i) is after 4>f{t) and (t>h{t) is before 4>b{t). Employing Eqs. (||) and (|l^), one 
may see that with this convention V^{t,t) = 'D^{t,t + St) = and V^lt^t) = 'D^{t + St,t) = 0. In the continuous 
notations it is equivalent to the choice of the regularization, where 9{0) = 0. As a result, for any number ^ of retarded 
or advanced functions one has 

Tr{X>^2?^...2?'^} = 0, (13) 

where the traces are understood as time integrations. If a problem is translationally invariant in time - the same 
statement has the simple reinterpretation 

oo 

J dujV^{Lj)v^{Lj + ni)...v^{Lj + n„) = o, (14) 



^For fermions it is more convenient to transform and ip in the different manner (it is allowed since ip and ip are not actually 
complex conjugated, but rather notations). As a result, the standard matrix form of the fermionic propagator is 

{ g^{t,t') J ■ 

In this notation fermion-boson interactions take an especially simple form Sint ~ 4'a'4''y°'^, with 7"^' = 1 and 7' = ai [ p^ . 
Since fermions are never a classical field, this form of the fermionic propagator is not related to our subsequent discussion. 

*This is true for a single function as well. For example, if T)^ = {Dq^ — iuj)~^ , one may worry that J duj{Dq^ — iio)~^ is 
divergent. This is not the case, however, since this integral is simply I)^{t,t) = 0. 
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which is an obvious consequence of retarded Green function being an analytic function in the entire upper half plane 
of complex energy. From Eq. (11) it is obvious that 



p^ = [2?^]^ (15) 

where Hermit conjugation includes transposition of the time arguments. 

In our simple example ~ —i{n + l)e~*'^°'^*~* and 2?^ = — me"*'^"*^*^* \ where n = n{ujQ) — p{ijJq)/{1 — p(u)f))) is 
the bosonic occupation number (since the system is non-interacting the initial distribution function does not evolve). 
As a result 

P«(i, t') = -ie{t - <')e-^""(*-*') ; 

V^{t, t') = -z(2n(cjo) + l)e-*"«(*-*') . (16) 

Notice that the retarded and advanced components are independent on the distribution function, whereas the Keldysh 
component does depend on it. This statement is much more general than our toy example. Going to the frequency 
representation, one finds that 'D^[lo) — —2TTi{2n{ojo) + l)S{oj — ojq) = —2Tri{2n{uj) + 1)S{uj — luo). In the particular 
case of the thermal equilibrium, one has 

V^{lo) = coth^ i'O^i^^) - T^^H) ■ (17) 

The last equation constitutes the statement |^ of the fluctuation-dissipation theorem (FDT) . It implies (as we shall see 
below) the rigid relation between response and correlation functions at equilibrium. FDT is always satisfied in the 
thermal equilibrium even in the presence of the interactions. We shall see later why the interactions do not change 
this relation. In general case it is convenient to define a function (matrix in the time space) F via 

and refer to F (or better its Wigner transform) as the distribution function. At equilibrium F{t — t') F{u!) — 
coth^ = 2n{uj) + 1, where n{uj) is the Planck distribution. 



6. Keldysh action and causality. Once we have understood the structure of the propagator in the cl ~ q space, we 
may go back to the action. Inverting the matrix Eq. (11), one finds 



A 1 

S[^,i,<j>,] ^2 J J dtdt' (Mt), Mt)) ( [p^] -1 ) ( t^^P) ) . (19) 

where [D~^]^ = [V^]-^ F ~ F[D^]-^ 7^ [V^']-\ Of course, in our simple example = S{t-t')[idt' -wo±«0] 

and [I?^^]-'^' — and therefore the correlator matrix is basically S{t — t')[idt' — t^^olci, as it should be. The point is 
that the structure of the Gaussian action given by Eq. (|l^) is absolutely generic and correctly encodes regularization 
of the functional integral. Since the Keldysh component carries information about the density matrix, there is no 
need any more to remember about the discreet representation. As we shall see below this structure remains intact 
under any perturbative renormalization. For the lack of a better terminology we shall refer to the robustness of this 
structure as the causality principle (since it is the reason for the response functions being casual, see below). Let 
us devote some time to this structure. 



^The analogous statement for fermions reads as 

g^(e)=tanh^ (0%) - g-\e)) . 

^Notice that the matrix in the exponent of the functional integral, iV~^ , is anti-Hermitian ! The convergence factors are 
hidden in the regularization. 
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Consider first the cl — cl element, which is zero. This zero may be traced back to the relation (^). It has, however, 
the much simpler interpretation. It says that if one takes a purely classical field configuration, 0^ = 0, the action is 
zero. But this is obvious, since in this case 0/ — (i>b and the action on the forward contour cancels exactly that on 
the backward. We therefore arrive to the simple but extremely important statement: 

-0] -0. (20) 

From this perspective the relation (^) is simply one of manifestations of the causality. Obviously Eq. (^o|) is not 
restricted to the Gaussian action. Consider e.g. the simplest interaction action Sint — ^ J^dt (0'/')^ = ^ Jdt [{4'f4'f)'^ ~ 
(^&0b)^]- Notice that this is the continuous notation; in discreet version (f> fields on the Keldysh contour C are taken 
at one step 6t later time with respect to (p fields. Transforming to the classical-quantum components, one finds 

oo 

S^nt [^cl 4A Sfty dt (</.2, + ^ (21) 

— OO 

which obviously satisfy Eq. (^0|). 



a) b) 

FIG. ^ First and second order 

diagrams for the {(fi(j>)^ theory. 

The fact that cl ^ q and q — cl components in Eq. (|l9| ) are advanced and retarded (means triangular in the time 
representation) is of the fundamental importance as well. Due to this fact Eq. (20) is intact in any order of the 
perturbation theory. Consider, for example, the first order correction in A to the cl — cl element of the Gaussian 
action (which is zero). The corresponding diagram is depicted on Fig. |a. From Eq. one reads out 



SSci-ci = 4A5Ry dt4>ci{t)V"it,t)ct),i{t) . (22) 

— oo 

As was explained above, 'D^{t, t) = 0, and the cd — cl element obviously remains zero. Employing Eqs. (p^), ( p^ ) one 
may show that all higher loop corrections do not change this fact. One may say that Eqs. (|l3|), ( p^ are part of the 
causality principle. To complete the prove that causality is preserved in perturbation theory, one has to show that 
cl — q and q — cl components remain purely advanced and retarded correspondingly. A reader is advised to check that 
this is indeed the case in the second order. Fig. ||b. 



7. Saddle point equations. Before doing the perturbation theory we have to discuss saddle points of the action. 
Eq. ( pO| ) states that there are no terms in the action which have zero power of both 0^ and ipq. The same is obviously 
true regarding 5S/54>ci and therefore the saddle point equation 

0(pcl 

may be always solved by 

$g = . (24) 
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The saddle point solutions are denoted by capital $. One may check that this is indeed the case for the action given 

ec( 

= = 0, (25) 



by Eqs. ([l9|), (21). Under the condition Eq. ( p4[ ) the second saddle point equation takes the form 

1 SS 



2 S(j)q 



where (^^[^c;] is retarded operator describing the classical dynamics of the field. In our example 

0''[<^ci]=tdt-ujo + X\^ci? (26) 



and the saddle point equation (25) is the time-dependent " Gross-Pitaevskii" equation. We have arrived to the 
conclusion that among possible saddle point solutions of the Keldysh action, there are always some with zero quantum 
component and the classical component obeying the classical equation of motion. We shall call such saddle points 
classical or casual (since perturbation theory on top of a classical solution preserves the causality structure). Notice, 
that the classical equation of motion is result of variation over the quantum component (in the limit where the later 
is zero). 

In view of Eqs. ( |20| ) and (|2^), the action on any classical solution is zero 0. Recall now that we do the theory 
for Z — 1 - this is the exact relation. To reproduce it in the saddle point approximation one has to show that the 
fluctuation determinant on top of the classical saddle point is unity. It may be a formidable task to demonstrate it 
explicitly, but it is helpful to remember that this must be the case. 

An extremely interesting issue is an existence and role of non-classical saddle point solutions, which have <I>^ ^ 0. 
Such solutions, if exist, have, in general, a non-zero action and may be loosely called instantons. Examples of 
instantons include (thermal) escape from a quantum well ||2^ and Wigner-Dyson statistics in metal grains . Since 
the normalization Z = \ '\s entirely due to classical saddle points, instantons must have strictly zero contribution to 
the partition function, Z . They may, however, contribute to correlation and response functions. The subject is poorly 
developed (unlike instantons in the imaginary time) and requires much more attention. 



8. Semiclassical approximation. Let me concentrate now on a classical saddle point, Eqs. (p4[), (|25|). To include 
quantum effects in the semiclassical approximation one has to include fluctuations of the 0, field near zero. To this 
end we keep terms up to the second order in in the action. The resulting semiclassical action takes the following 
general form 

oo 

dtdt' [^g[V-Y<f><] + {^gO''[cj),i]cf>,i + c.c.)] ; (27) 

c.c. stays for the complex conjugation. From this point one may proceed in two directions: 

(i) Since the action Eq. ( P7[ ) is Gaussian in 0g one may integrate it out and end up with the theory of the single 
filed, (pel- The corresponding action is 

oo 

Sscilcl^ci] = 2 / [ dtdt' M0^[^ci]V^)[V'']-\v''0''[cl>4)^,i . (28) 



If the non-linearity in (j)ci is neglected, one ends up with the Gaussian theory governed by the inverse Keldysh Green 
function, S ^ 4>ci['D^]~^4>ci f\ In general, the action Eq. ( p8| ) is pretty complicated. If it contains only a finite number 
of the degrees of freedom and is short-ranged in time, one may use the transfer-matrix technique to write down the 
corresponding " Schrodinger" equation. The result is the Focker-Planck equation. 



''This statement is a close relative of the Parisi-Efetov-Wegner theorem in supersymmetry or the fact that a replica symmetric 
saddle point results in zero action in the replica limit. 

8j£ 2,-R,A ^ (0 ^ ii^y^, then in the thermal equilibrium (c.f. Eq. (|l|)) , [D^]"^ = ^utanh^. In the T = limit one 
obtains the non-local quantum dissipative action, S ~ |ci;||<jic;(ii')P, while in the high temperature regime the action is local 
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(ii) One may perform the Hubbard-Stratonovich transformation with the auxiliary stochastic field ^ (i) to decouple 
the quadratic term in Eq. (|2^). The resulting action is linear in <j)q. Therefore the integration over (j>q leads to the 
functional i5-function. As a result one obtains stochastic Langevin equation ^ 

0''[UMt)^m, (29) 
where ^(t) is a Gaussian noise with the correlator 

mm) = '^[T^~Yit,t'). (30) 

For the diffusive dynamics, [P^'-^J^i = DV^ =F in the thermal equilibrium i['D-^]^{uj) = icoth^([P^]~i(w) - 
[I>^]-i(w)) 4T in the limit uj <^T. In this case one obtains the classical Langevin force with the correlator 
{(,{t)£,{t')) = 2TS(t — t'). The opposite limit is called sometimes the quantum Langevin equation. 

Notice that the above procedure may be reversed. Start from the classical Langevin equation, then introduce 
functional 5-function. Finally exponentiate it with the auxiliary field called 4)q and integrate out the Gaussian noise. 
Such procedure is known in the literature as Martin-Siggia-Rose-DeDominics (MSRD) Q technique. The result is the 
semi-classical Keldysh action, Eq. (|27|), in the high temperature limit. In the MSRD literature the field 4>ci is called 
(/), while (pq is usually called (j)- One may wonder if one can form their symmetric and antisymmetric combinations 
(j}f = (j) + (f> and 0b = — to recover an "underlying quantum theory" with the action S[4if] — S[(j)b]. In general 
it does not work since MSRD is only the semiclassical limit of Keldysh. Therefore terms of the order 0{(p^) may be 
" missing" . 



9. FDT and kinetic equation. Having understood the structure of the saddle points, we can construct a perturbative 
expansion around a classical saddle point 0. Let say we are about to calculate the full Green function, defined as 
Va^i3{t,t') — —2i{(j)a{t)(j)fj{t')), where a,(3 = cl,q. This is done in a standard way by expanding the exponent of the 
non-Gaussian action and applying the Wick theorem according to Eq. (|l^). One then rearranges the perturbation 
theory and defines the irreducible self-energy S through the relation: 

V^Vo+ Vo^Vo + Vo^VoVqY^Vq + . . . , (31) 



where multiplication is understood as the matrix one in the time and cl — q spaces; Vq is the bare propagator, Eq. (|l] 
The causality guarantees (this may be checked order by order) that the self-energy matrix has the structure of the 
correlator (i.e. V^^; cf. Eq. (|l|)) 

y (ff'\-f ^^{t,t')\ 



I]^(t,t') S^(t,t') 

Consider e.g. the second term on the r.h.s. of Eq. (^l|), after a little algebra one finds 

SVcl, cl = V^y^V^ + 2?^E^X>^ + V^T.'^V^ ; 
5Vd, q = V^T.^V^ ; 
5Vq, cl ^ V^y^V^ ; 

SVq,q = 0. (33) 

As a result the q — q component remains zero, while cl — q and q — cl remain retarded and advanced correspondingly. 
Indeed, product of triangular matrices is again a triangular matrix. Obviously, the same is true in all higher orders as 
well. Therefore the causality structure remains intact even in the presence of non-linearity (at least in perturbation 
theory). Eq. (p3l|) may be rewritten as the matrix Dyson equation 



®One may worry about appearance of the Jacobean upon integration of the functional 5-function over It does not happen 
with our regularization of the functional integral (Jacobean is unity). The reason is exactly the same which guarantees Z = 1 
upon integration of the fluctuations around a classical saddle point. 

^''I am not aware about any attempt of doing the perturbation theory on top of a non-classical saddle point. 
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- E) P = 1 



(34) 



It is a simple matter to show that the q — cl component of this equation is nothing but the kinetic equation on the 
distribution function F. Since this procedure is well documented in the literature [0j^J^ we shall not describe it 
here Let mc only mention that to obtain the textbook (semi-classical) kinetic equation one has to assume that 
all external fields are slow on the scale of temperature. This assumption justifies the time-Wigner transform of the 
matrix F, which is a necessary step in the derivation of the kinetic equation. 

The other way to obtain the kinetic equation is to introduce time non-local field Q{t, t'), canonically conjugated to 
the composite operator (f)(t)(j){t'). A small amount of disorder makes Q{t,t') a dynamic field which is governed by a 
certain effective action (this action may be obtained by integrating out the fields (j) a-nd 0) . The saddle point equation 
for this effective action Sef / [Q] is the kinetic equation. As usually, the (semi) classical (kinetic) equation is obtained 
by variation over the quantum component of the Q-field. For the details of this procedure see Ref. One can 

possibly go beyond the saddle-point approximation by including the Gaussian fluctuations of the quantum component 
of Q (pretty much in the way we did above) . The result of such procedure is probably the Langevin-kinetic equation 
of Kogan and Shulman ||2|]. To the best of my knowledge, it was not demonstrated explicitly yet. 



10. Sources, external fields, response and correlation functions. Up until now we have been busy with complicated 
representation of unity: Z ^ I. This is an interesting, but useless job. To make the technique work we need to 
calculate something more informative. This something is the generating function, which is obtained by adding the 
source fields to the action. It is convenient to discuss it in parallel with introduction of external classical fields (such 
as scalar or vector potentials). Let us add the source term Ssource = \ J^dtVcfxj) to the action. The source field V{t) is 
defined on the Keldysh contour, therefore it is rather a couple of fields (V/, Vb), residing on the forward and backward 
branches correspondingly. We decompose it now on the classical and quantum components exactly in the same way 
we did it with the (j> field: Vd = {Vf + Vb)/2 and Vq — (V/ -- Vf,)/2. The source action takes the form 



'^^The simplest exercise is the theory of the real bosonic field with the interaction Xcj>'^ — > 3X4>ci4'ci + In the second order 
the self-energy Fig. | is given by E^-^ - X^V^ * Vq'"^ and T,^' ~ A^(r>^ * + * 2?^ + * V^), where the star denotes 
convolution in the energy space. 




a) b) 

FIG. H Example of the Keldysh diagrammatic for the tj)^ theory. Fields cj>ci are denoted by full lines and 0, by dashed lines. 
There are two vertices with relative coefficient 3. There is one retarded, one advanced and three Keldysh self-energy diagrams 



Employing the "magic" identity cothecoth(e + uj) — 1 = cotho; (cothe — coth(e + uj)), one may immediately see that if Do 
satisfy FDT, Eq. (0, then so does E, namely T,'^ = coth^(E^ - E^). Then the Dyson equation (||) says that the exact 
propagator also satisfy FDT. The "magic" identity is directly responsible for the fact that the Planck distribution (coth^) 
nullifies the collision integral. 
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oo 




— oo 



Notice that the quantum source is coupled to the classical component of the density (defined in a usual way as a half 
sum of the density on the forward and backward branches) and vise-versa. We define now the generating function as 

Z[Vci,Vg]^{e'''—) . (36) 

The classical component ol the source field is nothing but classical external scalar potential, which is supposed to be 
the same on the two branches. The presence of such field does not violate the fact that the forward and the backward 
evolution brings the system back to the initial state. Therefore it does not change the normalization Z = 1. We 
therefore arrive to the conclusion 



Z[V,i,Vg = Q] = 1. (37) 

One can also check this statement explicitly by e.g. expanding the action in powers of Vd and making contractions. 
All terms in this expansion vanish due to the fact that {<^q4>q) = 0. Eq. ( |37| ) is a close relative of Eq. ( |20| ) and is just 
another manifestation of the causality. As a result, adding to the theory classical fields alone does not do any good - 
we still work with unity. Therefore, one necessarily has to introduce quantum sources (they change sign by going from 
the forward to the backward branches of the contour). The presence of such source fields explicitly violates causality, 
and thus changes the partition (generating) function. On the other hand such fields usually do not have a physical 
meaning (see, however, the next section) and play only auxiliary role. In most cases one uses them only to generate 
observables (by appropriate differentiation) and then put them to zero, restoring the causality structure of the action. 
Notice that the classical component, Vd, does not have to be taken to zero. 

Let us see how it works. Suppose we are interested in the average density of bosons at time t in presence of a certain 
scalar potential Vci{t). Of course, we are looking for the classical component of the density. According to Eqs. (|35| ) 
and (|3^ ) it is given by 

n{t; Vd) = -I TTTTT^ZlVd, Vq] . (38) 



svsr 



v„=o 



To calculate this quantity one has to find first an appropriate (classical) saddle point of the action in presence of Vd 
(as was mentioned above this amounts to solving the kinetic equation) and then study quantum corrections. The 
problem is simplified if the external field, Vd, is weak in some sense. Then one may go away with the linear response 
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2 



5Vd{t') ' ' ' K,=o 5Vd{t')5Vq{t) 



I put the subscript R, anticipating on the physical ground that the response function must be retarded. We shall 
demonstrate it momentarily. First, let us introduce for esthetical reasons the matrix 

The fact that Xci.ci = is obvious from Eq. (^7|). To evaluate the linear response matrix x consider first the Gaussian 
action, Eq. (p^. Adding the source term, Eq. (|35|), and integrating out the boson fields 0, one obtains 

Z[l/] = e-Tnn(i+w) ^ (41) 

where V — [VqCfQ + VdUi)/'^- Notice that the normalization is just right, since Z[0\ — 1. One may now expand 
ln(l+2?F) to the second order in V and then differentiate twice. It is a straightforward matter to see that the causality 
structure of 2? imposes it also on x- E.g. X^it,^') ^ TT{VaiVao} = V^{t,t')V^{t',t) + {t,t')V^{t' ,t), which is 
obviously retarded. In presence of interactions, Eq. (pT]), one may perform Hubbard-Stratonovich transformation and 
then integrate out 0, </> fields. At this stage the causality structure of x becomes apparent (it does not mean yet that 
one knows to evaluate x)- 

The purpose of the above discussion was to demonstrate that even for the linear response problems (which may be 
always solved in Matsubara technique) the Keldysh formulation is natural and transparent. For one thing it allows 
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one to overcome the analytical continuation procedure, which may happen to be a headache. The other advantage is 
that it gives both response and correlation functions in the single framework |^. 



11. Counting statistics and counting fields. A very interesting application of the Keldysh formalism was suggested 
recently by L. S. Levitov and coworkers Suppose we are interested in a certain quantum mechanical observable 
O (for example, a total charge which passed through a wire during a long time tf ~ ti). One may look for e.g. its 
(quantum) expectation value, (O), or noise (O^) — ((O))^, or even higher order moments. The most general quantity, 
however, is the distribution function: P{0) - probability to observe the value O in a single measurement of O. If one 
knows P{0), then (O*^) = / dOP{0)0^ . Instead of discussing P{0) it is convenient to consider its Fourier transform 
F(A) = / dOP{Q)e^^^ . The variable A - is called the counting field, and P{0) (or F{\)) is called the counting 
statistics. 

From the above expressions one may deduce that 

' i\6 



F{\) ^ i^e'^^j . (42) 

The non-trivial fact is that the angular brackets in this expression must be understood as the averaging along the 
Keldysh contour. The prove of this statement requires a careful consideration of the time-ordering convention for the 
higher moments of O. We refer a reader to the orig inal works 0jl| for more details. If the Keldysh nature of the 
averaging in Eq. ( ^^ is accepted, it should be obvious by now that the counting field A is a purely quantum field. 
Indeed, A was introduced to generate moments of O upon differentiation. As we have seen in the previous section - 
it is the quantum component of the source field (having different sign on forward and backward branches (!)), which 
serves for this purpose. As a result, Eq. ( ^ ) suggests that the Keldysh action of the system at hand must be modified 
as 



S^[^]^S[^]+ J dt\{t)6ci{t). (43) 

— oo 

We put O — j^' dt Ociit), where the classical observable dci(t) is expressed through the fields 4>{t). We have also put 
\q{t) = X9{t — ti)9(tf — t) - thus the name - counting field. 

The novel feature is that the quantum source, Ag, should not be taken to zero in the end of calculations. It rather 
remains finite and possibly even large to allow the inverse Fourier transform of F{\). The presence of the finite 
quantum source in the action spoils immediately the causality structure of the theory. For example, S'^[4)q = 0] ^ 
A ^ 0. As a result, <I>g = is not a saddle point solution any more. There are no classical saddle points in the theory 
! A new non-classical saddle point has in general a non-zero action Ssp = S'sp(A) ^ 0. Neglecting fluctuation effects, 
one finds for the generating function of the irreducible moments: —i\n.F{\) — Ssp{X). Such program (including 
fluctuations) was realized by Yu. V. Nazarov |lj] for the counting statistics of charge passing through a diffusive wire. 
See also related works on counting statistics of adiabatic pumps p^ . 

To conclude this section I want to mention a striking similarity between counting statistics in Keldysh theory and 
statistics of rare events in supersymmetric technique p3|. To my opinion this subject is far from being closed. 



12. Keldysh technique in disordered systems. Due to the presence of the intrinsic normalization condition, Z = 1, 
the Keldysh technique is ideally suited for treatment of disordered systems, such as spin glasses or dirty metals. 
Notice that the response and correlation functions were defined by differentiation of Z[V], Eq. (|40|), rather than more 
familiar In Z[V] (since Z[Vq = 0] = 1 this is completely equivalent). Let me also mention again that the normalization 
constant, Tr{po}, introduced in Eq. (^), is disorder-independent. The point is that disorder is switched on (and off) 
on forward (backward) parts of the contour after (before) t — —oo. As a result, the "partition" or generating function 
Z (rather than its logarithm) is the proper object to be averaged over the quenched disorder. This circumvent the 



^^In the equilibrium response and correlation are connected by the FDT; i'^) ~ coth^(x^(tt^) — x^i'^))- However, out of 
the equilibrium they require separate calculations. 
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need to introduce replica (or supersymmetry |25|) to take care of the logarithm. This was first realized in the 
context of classical spin glasses (MSRD version of Keldysh) by H. Sompolinsky |Q and substantially advanced later 
. In the context of electrons in the random potential the corresponding theory (Keldysh non-linear cr-model) was 
developed in Refs. |15 



Here I restrict myself to a few remarks only. On the perturbative level the replica and supersymmetry methods are 
supposed to eliminate "parasitic" diagrams. Such decoupled closed loop diagrams are generated if one erroneously 
averages Z instead of InZ. The replica trick first prescribes factor n (number of replica) to such loops and then 
eliminates them by taking the limit n — > 0. In the supersymmetry bosonic loops exactly cancel fermionic ones. The 
Keldysh analog of such cancellation mechanism is given by Eqs. (^ and (p^. Indeed, the closed loops coupled to 
the rest of a diagram by the static disorder carry either only retarded or only advanced Green functions |^ The later 
vanish according to Eqs. (^3|) (p^). Therefore it is not surprising that all perturbative results of the non-linear cr-model 
(including electron-electron interactions) are correctly and elegantly reproduced in the Keldysh version p5| , p6[ |. 

Whether or not the Keldysh cr-model contains the non-perturbative effects (such as instantons) - is a hundred 
dollar question. I believe that the answer is positive. At least, it is possible to obtain the non-perturbative (though 
only asymptotic) part of the level statistics via the Keldysh cr-model |l^. To this end, one has to analyze carefully 
non-classical (having a non-zero action) saddle points on the massless manifold of the cr-model. In any case, a lot to 
be done before the Keldysh technique will become a reliable tool for the non-perturbative calculations. 



II. DOI-PELITI TECHNIQUE FOR CLASSICAL SYSTEMS 

1. A bit of history. The technique was invented by Doi in 1976 and later rediscovered by Peliti ||T^, and others 
p7| , p8[ . J. Cardy and B. P. Lee were probably most persistent in applying the technique to various non-equilibrium 
statistical problems p9|. These notes are partially based on lectures of J. Cardy pO(|. 



2. Motivation. We are talking about extremely broad and diverse types of problems whose applicability ranges 
from epidemiology and dynamics of biological populations to models of chemical reactions and glasses. The problems 
at hand are usually determined on some d-dimensional lattice. They are defined by specifying the following entities: 

• The agents, usually named as A, B , . . . These may be molecules, animal species, stock options, vehicles, etc. 

• Rules of agent's motion on the lattice. It may be e.g. random walk, or directed ballistic motion. 

• Rules of agent's interactions. E.g. the simplest one is the binary annihilation: if two agents meet at the same 
lattice site they annihilate with probability /i. The short way to write is 

A + A-^0. (44) 

This model may, for example, describe a simple chemical reaction, where two identical molecules form a stable 
aggregate which segregates from the solution. An other famous example is the Lotka-Volterra system [ pT|p2[ 



B^2B- (45) 



A + B^2A, 



where A designate "foxes" and B — "rabbits" . Foxes die with probability /i; rabbits proliferate with probability 
cr, if a fox and a rabbit meet — the rabbit is eaten and another fox is produced with probability A. 



^^Indeed, the static disorder Udis enters the Keldysh matrix action as Udis<^i- Integrating out of the fields one ends up with 
Tr ln(l + VaiUdis), cf. Eq. (kll). Expanding the logarithm one finds terms of the form TrjOcri . . . Vai}, the later are obviously 
TrjC^ . . . V^} + Tr{V-^ . . . O*} = 0. 
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3. Mean-field. The simplest way to treat systems described by Eqs. ( |44[ ) or ( [45| ) is to employ the mean field 
equations. To this end one introduces the mean density of agents, say A, at a lattice site r at time t - n^(r, t). For 
the system Eq. (^4|) with random walk motion (ARW - annihilating random walk) the mean field equation for the 
density is 

^^DVluA -t.nl, (46) 

where D is the diffusion rate. Derivation of the mean field equation is self-evident from the definition of the model, 
Eq. (0). We shall see later, that from the perspective of Doi-Peliti technique the mean-field is nothing but equation 
for a classical saddle point. A reader is advised to write down and investigate the mean-field equations for the 
Lotka-Volterra system, Eq. (|4|) . 

The mean-field ignores fiuctuations. It works reasonably well, when the number of agents at each site is large. 
In the opposite limit where the average number of agents on a site is less then unity, the system's behavior may 
be solely determined by fluctuation effects (especially in low spatial dimensionality). For example, Eq. (|4^) predicts 
UA ^ t^^ at i OO, while the correct result ]2^ , |30t is nA ^ t~'^/'^ for d < 2 and ua ^ t^^ for d > 2. The behavior is 
clearly very different from the mean-field prediction for dimensionality below the critical one, dc = 2. A traditional 
approximate way to handle fluctuations (which also stops short to obtain the above mentioned results) is by means of 
the Focker-Planck or Langevin equations. I shall comment below on the approximations involved in their derivation. 



4. Master equation. The formally exact way to handle the problem is the master equation. One first defines a 
many-body microscopic configuration a, determined by the (integer) number of agents at each lattice cite 

a = {ni,n2, . . . riAf} , (47) 

where M = i'' is the total number of cites in the d-dimensional lattice. Next one defines P{a, t) - probability to find 
the system in the many-body configuration a at time t. The master equation has the form 

^F(a, t)^Yl T^P^c^PiP, t) " ^ 7^„^^P(a, t) , (48) 

where the sums run over all many-body configurations /3. The first term on the r.h.s. describes "in" , while the second 
one "out" processes. The coefficients TZp^a and TZa^p arc the transition rates. A particular case when, 

7^a^/3 Wia) ' ^^^^ 

where Vl^ is a certain positive function uniquely defined for each state, is said to satisfy the detailed balance condition. 

As an example let us consider the binary annihilation on a single site. The microscopic state is specified by a 
positive integer number of agents, a — {n}. According to Eq. ( |4^ ) the master equation takes the form 

„/ N (n + 2)(n + 1) , nhi — 1) , , , 

-Pin, t) = fi ^ >j Pin + 2,t)-ti ^ ^ > Pin, t) . (50) 

The only non-zero rates are 7?,„+2-»n and which are proportional to the number of possible pairs in the 

corresponding state. These rates do not satisfy the detailed balance. 



5. Quantization. The basic idea Q is to draw the analogy between the master equation and the many-body 
Schrodinger equation. Indeed, both equations are linear in P (or ^E*) and first order in df. To make the analogy 
explicit one introduces the "second quantization". Define the ket-vector \n) as the microscopic state with n-agents. 
Let us also define vector 

oo 

!*(<)) ^^P(n,i) In). (51) 

Note that the weight, P, is probability rather than the amplitude. In these notations the master equation takes the 
form 
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I I* W) = E = f £ + + l)P(n + 2, t) - nin - 1)P(,^, i)] 1") ■ (52) 

n=0 n=0 



The aim now is to write the r.h.s. of this expression as — where H is the "Hamiltonian" operator. To this 
end we introduce the creation and annihilation operators: 

aV) = l"+l); (53) 
a\n) = n\n - 1) . (54) 

As a byproduct, one has a|0) = 0. One may immediately check that such operators are "bosonic" ^ 

[a,a^] = l. (55) 

As for any pair of operators satisfying Eq. (|55| ) one may prove the identity 

e°/(a,at) = /(a,at + l)e", (56) 

where / is an arbitrary operator-value function. With these definitions one may check that the Hamiltonian, we are 
looking for, has the form 

i^=f((«^)V-a2), (57) 

where the first term on the r.h.s. is the "out" and the second one is "in" term. Generalization for the many site 
problem with the random walk (ARW) is straightforward, the resulting Hamiltonian is 

H = Dj2ial- a])ia, a,) + | ^ " ' (58) 

<ij> i 

where < ij > are nearest neighbors. 

6. Observables. So far we succeed to put the master equation into the Schrodinger form 

l\M,it)}^-H\^it)}. (59) 

The differences with the true Schrodinger equation are: 

• There is no imaginary unity, meaning that we are dealing with the imaginary time Schrodinger equation. 

• H is, in general, non-Hermitian (unless the detailed balance condition is satisfied). 

• The most important: the mean value of some observable O, is not {'^{t)\0\'^{t)) . Indeed, such combination is 
quadratic in P, while the mean value should be linear in probability, P. 

To obtain the proper mean value one defines the bra-coherent state 

(^'o| = (0|e^ (60) 



This state has the special property: (4'o|'i) = 1 for any n. This fact may be checked using Eq. (p6[). Employing this 
coherent state, one may show that the mean value of an observable may be written as 

{0) = {^o\d\^it)) . (61) 
Indeed, let us consider density as an example, O — a'' a 



In problems where not more than one agent is allowed on a given site, (e.g. traffic flow models) one may probably employ 
"fermions". I am not aware, however, of such attempts. 
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{n{t)) = {^oWa\^{t)) = (*o|a|*(t)) = (0|eVa ^ P(n, t)!") = , 

n n 



(62) 



in agreement with the expectations based on the meaning of P{n,t). The second equahty on the r.h.s. is very 
important; it is an immediate consequence of the coherent state being the left eigenstate of the creation operator, 



7. Normalization. To have a consistent theory, one has to worry that the total probability to find the system in 
some state is unity 

l = ^P(n,0-(*o|*(t)). (63) 

n 

The last equality is just the mean value of the unit operator. One may solve formally the Schrodinger equation 
and write |*(<)} = exp{-iJ(at, a)t}|^'(0)). Eq. (|6|) is satisfied iff {^n\H{a\a) = 0. Since the coherent state is an 
eigenstate of the creation operator, (^Pola^ = (^o|, one arrives at the conclusion that any legitimate Hamiltonian must 
obey 

iJ(at l,a) = 0. (64) 

E.g. the Hamiltonian of ARW, Eq. (^8|), indeed satisfy this condition. To introduce the functional integral represen- 
tation we'll need to have the normal ordering of creation and annihilation operators. Therefore it is a good idea to 
commute e° through the Hamiltonian employing Eq. (pq). 



(0|e°e-^('^''")*|^'(0)) = (0|e--^('^*+i^'^)*e"|*(0)) . (65) 

If an observable is considered, the same commutation may be performed provided that the observable is written in 
terms of annihilation operators only. An example, of this procedure is provided by the second equality in Eq. (|6^). 
We shall redefine from now on the Hamiltonian as H{a\a) = H{a^ + l,a). In terms of the new Hamiltonian the 
fundamental normalization condition, Eq. (|6^), is satisfied iff 

H{a) ^Q,a)=Q. (66) 

We shall see momentarily that this condition is the direct analog of the Kcldysh normalization condition, Eq. (pO|). 



Field theory. Consider the "partition function" of the theory 

Z= (0|e--f^('^^'^'^)*|*(0)) 



(67) 



where |5'(0)} — e" |5'(0)} is the renormalized initial microscopic state, its precise form is not very important for our 
further discussion. Due to normalization condition and according to Eqs. (|63|), (|65|), Z = 1. We divide now the time 
interval [0,t] onto N ^ oo slices and introduce the coherent state resolution of unity at each time slice. As a result, 
one obtains in the standard way [|l9| 



Z= /P0*0exp{-5[0*(O,0(O]}, 



where 

For example, in the case of ARW one has 



dt 



Hi 



(68) 



(69) 



S[4>*{t,r),ct>{t,r)]^ jdt ldr[r{dt-DVl)^+^^{iq^*)^^'+2q^*4>^) 



(70) 
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9. Analogy with the Keldysh technique. Notice that both in Keldysh and Doi-Peliti techniques we deal with the 
"partition function" Z = 1. (No matter that the reasons for this normalization are very different: closed time contour 
in Keldysh versus conservation of probability in Doi-Peliti.) In both cases we succeed to write the partition function 
as the coherent state functional integral, cf. Eq. (|^) and Eq. (^) . To this end in both cases we had to double number 
of degrees of freedom: classical and quantum in Keldysh versus (j) and 0* in Doi-Peliti |^. The normalization requires 
the strong condition on an acceptable action: Eq. (Eoh in Keldysh versus 



5[r=O,0]=O (71) 



in Doi-Peliti. This equation is a direct consequence of Eqs. (|66|) and (|69D. Comparing Eqs. (|2C|) and ( [7lD one may 
suspect the analogy 

0* ~ 09 ; ~ ■ (72) 

We shall see below that this analogy indeed works for all intends and purposes. Let us rewrite the quadratic part of 
the action of the ARW model in the following form 



-dt - DVl 
dt - DWl X 

y — ^ V / 



(73) 



Here R and A designate retarded and advanced linear classical operators. Symbol K stays for the Keldysh component 
of the correlator. In our simplistic example it is absent, X = 0. This fact is not generic by any means: in more 
complicated models e.g. Lotka-Volterra there is non-zero Keldysh component. Moreover, even if it is zero in the bare 
action - it is usually generated in the process of renormalization. On the other hand, the fact that the — element 
is zero is not accidental. This zero is direct a consequence of Eq. ( |7l| ) and it is protected against any renormalization. 

Once the structure of the quadratic action is understood, one may write down propagators to build the perturbation 
theory (cf. Eq. (|ll]) ) 

{4>{tW{t')) ^V^itX) ^ {Dq^ - lu;)-^ 
{^*{t)<P{f)) ^V^{t,t') ^ {Dq^ ■ 
{^*it)<P*{t'))=0; 

(0(t)0(t'))=2?^(t,O- (74) 

As in the Keldysh case, 'D^{t, t) = T>^(t, t) = 0. It is exactly the same causality structure, which is responsible for the 
conservation of probability (normalization) and the entire internal consistency of the theory. The analogy between 
the two techniques goes actually much further than the perturbative expansion. 



11. Saddle point equations. Exactly as in the case of Keldysh, the saddle point equation 
may be always solved by 

$*=0. (76) 



The saddle point solutions are denoted by capital $. Under the condition Eq. ([76|) the second saddle point equation 
takes the form 



Notice that since we are interested in real density of agents, one may naively expect to obtain a field theory of the single 
real bosonic field. We have to deal with both and 0*, however, despite of the fact that the object is purely classic and the 
phase has no evident meaning. 
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(77) 



where 0-'^[(f>] is the retarded operator. In the example of ARW, Eq. ( ^0| ) 

O'^m dt - DWl + fi<i> (78) 

and the saddle point equation ( fn\ ) is nothing but the mean-field equation, Eq. (p6|). Notice, that the classical (mean- 
field) equation of motion is a result of variation over the "quantum" (</)*) component (in the limit where the later is 
zero). In view of Eqs. (^Tj) and ([76|), the action on any classical solution is zero. 

The existence and role of "non-classical" saddle points (with $* 7^ and S^p ^ 0) in the context of Doi-Peliti 
technique, to the best of my knowledge, was not yet investigated. 



12. Semiclassical approximation. To include fluctuation effects in the semiclassical approximation one keeps terms 
up to the second order in (j>* in the action. For the ARW problem the entire action is semiclassical (this is not generic 
by any means). Let us rewrite it in the following generic semiclassical form 



dt 



'[V 



-1-\K, 



(79) 



For the ARW problem [D^^]^ — ii(f)^ > 0. Positive definiteness of this quantity is necessary for the convergence of 
the (f)* integral. From this point one may again proceed in two directions: 

(i) Since the action Eq. (79) is Gaussian in 0* one may integrate it out and end up with the theory of the single 
filed, (j). The corresponding action is practically identical to Eq. (|2^). Employing the transfer-matrix method, one 
arrives at the Focker-Planck equation. 

(ii) One may perform the Hubbard-Stratonovich transformation with the auxiliary stochastic field ^ (t) to decouple 
the term quadratic in (jj* . Subsequent integration over </)* results in the functional 5-function. This way one obtains 
stochastic Langevin equation 



where ^(t) is Gaussian noise with the correlator 



(80) 



(81) 



Some studies start from postulating Langevin (or Focker-Planck) dynamics for problems at hand. We note that, 
this approach is not exact, but rather a semiclassical approximation. In some cases, like the ARW, the semiclassical 
approximation and hence the Langevin dynamics is indeed exact. However, even in such cases the correlator of noise 
is far from being trivial. Moreover, the noise term in Langevin equation, Eq. (|80|), contrary to any intuition, happens 
to be pure imaginary. It is hard to imagine how such dynamics could be postulated from the outset. 



13. Sources. To calculate an observable, e.g. mean density, one need to introduce sources. As was mentioned 
above, mean value of any observable may be expressed through annihilation operators only, c.f. the second equality 
in Eq. (^). In the functional integral language it means that pre-exponential factors at hand are given in terms of 
the "classical" field, </>, only. Therefore the source term, one needs, has the form / dtV* {t)(j){t) and the generating 
functional is given by 



Z\V 



[V*] = (e^ dtv (m(t)\^ 



(82) 



where the angular brackets stay for averaging with the action, Eq. (|69|). With this generating functional one may 
obtain, e.g., the mean density of agents 



it)) = (*o|a|*(t)) = /l?(/)*(/)e-^[*'(*)'"^(*)l(/)(0 = 



5 

5V*{t) 



Z[V* 



V'=0 



(83) 
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Other observables are evaluated in the similar fashion. Note that as in the Keldysh technique, observables are given 
in terms of the "classical" field, (/), while the generating function is functional of the fictions "quantum" source, V*{t), 
which is to be taken to zero after differentiation. 



14. Short summary. The aim of this chapter is not to give a detailed account of a specific calculation employing 
the Doi-Peliti technique. Lectures and papers of J. Cardy et al. [^,|9j are highly recommended for this purpose. 
My intention was to highlight the similarities and common concepts with the Keldysh techniques. The hope is to 
facilitate intermixing of the two largely non-intercepting communities. Let me stress once again common points of 
the two techniques: 

• Internal fundamental normalization, Z ^ 1. 

• Doubling of degrees of freedom: classical- quantunri!\ 

• Causality: retarded-advanced structure of the action and the propagators. 

• Existence of the classical extrema: = ; O^^ci — 0. 

• Langevin and Focker-Planck approaches as the semiclassical approximation — second order in (j)q. 

• Generating function of classical observables is functional of the quantum source. 

I am indebted to V. Lebedev, A. Andreev, A. Altland and V. Elgart for teaching me various aspects of the techniques 
and to L. V. Keldysh, M. Doi and L. Peliti for inventing them. This research was supported in part by the BSF grant 
N 9800338. 
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